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An orbital ordering occurs in many transition metal compounds with Jahn- Teller ions (Cu^"*", 

. . . Mn^^, low-spin Ni^"*", Ti'^"'' etc.). It plays an important role in these materials. At the same time, 

' exchange interactions in orbitally degenerate systems are inherently frustrated, even in materials 

' with simple crystal lattices. We discuss the origin of this frustration, considering in detail materials 

, with a cubic and triangular lattice of transition metal ions. We also discuss possible types of ground 

. states of such systems, e.g., disordered orbital liquids and ordering due to the order- from-disorder 

$^ ' mechanism. 
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I. INTRODUCTION 



Many transition metal (TM) compounds show an orbital orderingiiSi^. Typical examples are TM oxides with so- 
^ . called strong Jahn- Teller (JT) ions, e.g. Cu^+(d^), Mn'^+(d^), low-spin Ni'^+(d^) etc, which in the cubic crystal field 
^' ^ ' of a regular oxygen octahedron have one electron or one hole on the doubly degenerate Cg level. According to the 
well-known Jahn- Teller theorem, such states are inherently unstable, and at low temperatures a distortion of the 
local coordination lifts this orbital degeneracy. In concentrated systems it leads to a phase transition - the so-called 
' cooperative Jahn- Teller effect or orbital ordering^. An orbital ordering is also observed in materials with partially 
filled three- fold degenerate t23-levels of TM ions, e.g. Ti^+((i^) or V^+((i^). Some of these materials have very unusual 
propertiea^i^. However, the JT effect for t2g-electrons is usually weaker than for Cg-electrons and can be strongly 
,: affected by the spin-orbital interaction AL • S, which for Cg-electrons is essentially quenched. This complicates the 
Ch discussion of orbital ordering in t2g-systems and in this paper we concentrate on materials with eg electrons, though 
O many of our conclusions can be applied to t2g-systems too. 

An orbital ordering can be described similarly to a magnetic (spin) ordering. In the case of one electron or hole on 
a doubly-degenerate orbital, one can describe the orbital occupation introducing a pseudospin (or isospin) T = i, so 
that, e.g., the d^^^i _^'2 / d^'^ -y2 orbital corresponds to the state with = +^/T^ — A superposition of these two 
orbital states 
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O. \e) = cos-\3z'-r')+sm-\x^-y'), (1) 
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i: is associated with a particular point on a circle in the (T^, T'^)-plane, see FigQ] We only consider here linear 

' combinations with real coefficients. In principle, one can also consider states with complex coefficients, e.g., 

|: \±) = ^{\3z'-r')±^\x'-y')), 

. but they are rather special: they do not couple to the lattice, but break the time-reversal symmetry and have a 

' nonzero magnetic octupole momenlsi. Such states are eigenstates of the operator that usually does not appear in 

O the exchange Hamiltonian of doubly-degenerate JT systems. 

O The pseudospin operators can be used to describe interactions between orbitals on different TM sites in the same 

^ , way one describes the spin exchange. An orbital ordering corresponds then to a nonvanishing average pseudospin (T^). 
However, the orbital exchange usually has a more complicated form than the spin exchange, due to the directional 
nature of orbitals. The strong anisotropy gives rise to a frustration of orbital ordering even in materials with simple 

: bipartite lattices. 



II. PHONON-MEDIATED ORBITAL EXCHANGE 



We illustrate this point by considering effective orbital interactions in a perovskite lattice - a simple cubic lattice 
of TM ions with oxygens located on the edges between them (we ignore here for simplicity the tilting of (TM)06- 
octahedra present in many materials, which lowers the symmetry from cubic to orthorombic). There are two mecha- 
nisms of orbital ordering: a lattice instability due to the JT coupling of electrons to lattic o^'^'^ , and a purely electronic 
mechanism resulting from the virtual hopping of TM electrons between neighboring sitesi^. Both lead to essentially 
the same type of effective exchange interactions. We first discuss the lattice-mediated orbital exchange. 



FIG. 1: The {T^,T^) pseuodospin plane. The points A, B, C and D on the circle correspond, respectively, to the d:j^2_r2, 
d^2_y2, dg^2_^2, and d^y2_^2 orbitals. 
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FIG. 2: A pair of TM ions (black circles) along the z axis in the perovskite lattice. Arrows indicate the shifts of oxygen ions 
(white circles) resulting in lifting of the double degeneracy present for a regular Og-octahedron. 



Consider a pair of neighboring TM sites j and j + z along the z axis (see Fig|51l. For this pair the electron-lattice 
interaction responsible for lifting the double degeneracy of the e^-states has the form 

KQl 



HjT = -9 J2 {4lAl.-42A2^)Q3^ + ^, (2) 



where dli^ creates the state jSz^ — r^) at the site i with the spin projection a and dlj^. is the creation operator of 
the other e^-state — y^). The Q^i phonon coordinate is the linear combination of the shifts of oxygens (from the 
Oe-octahedron containing the TM site j) along the corresponding TM-TM bonds (see Fig[2J): 

" - '^j-z/2) - {uj+x/2 - - {uj+v/2 - ■Uj-y/2)] ■ 

In the subspace of states with one electron per TM site the combination of the electron operators and diaa in 
Eq.(|2Jl is equivalent to the pseudospin operator: 

J2[dl^da.-4,,d,2.)^2Tt. (3) 

(T 

Excluding the coordinate Uj^z/2 of the oxygen located between the two TM ions, we obtain for an effective orbital 
interaction for the pair 

H,,+. = JT^T^+., J=%, (4) 

which has the (antiferroorbital) Ising form favoring the antiparallel alignment of neighboring pseudospins, e.g. Tl = 
+ j (the d'iz'^-r^ orbital on the site j) and 7?^^ = — i (the dx2^y2 orbital on site j + z). 

However, complications arise when we consider effective orbital interactions in the whole cubic lattice of TM ions. 
Since in the cubic lattice the three directions x, y, and z are equivalent, we can describe the exchange along the x 
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direction by replacing the basis of states (|3z^ — r^), \x'^ — 'ip)) by (|3x^ — r^), — z^)) and similarly for the y direction. 
The symmetry of the cubic lattice under circular permutations of the x, y, and z axes is reflected in the three-fold 
symmetry of the (T^, r^)-plane (see FigQ, where the states d^^2_^2, d3x2_r2, and dzy2_^2 correspond, respectively, 
to 6* = (point A in Yig^ , 6* = +^ (point C), and 6* = - ^ (point D), where 

\e^±'-l) = -\\2.z^-r^)±^\x^-f). 
One can also introduce the corresponding symmetric notation for pseudospin operators 

so that — r^) = -|-^|37^ — r^), where 7 = x,y,z. Using the symmetry of the cubic lattice and Eq.Q, we can 

then immediately write the Hamiltonian of the exchange along the x axis. 



and similarly for a pair along the y axis: 



= Jl^r^^y = J y-\T: + ^T^j l^-^TJ + ^T/j . (7) 
The full Hamiltonian of lattice mediated effective orbital interactions has the form 

i i 

It includes also the on-site term (the second sum in Eq.||SJl), which suppresses complex orbitals discussed above by 
favoring alignments of pseudospins in the (T^,T^) plane. 

We can now see the source of the frustration in the orbital sector: while the exchange along the z direction stabilizes 
alternating d^z^-^^ and d,j.2_y2 orbitals, the equivalent couplings along the x and y directions favor other pairs of 
orbitals, i.e. [d^r^2_^2,dy2_z2 and {d^y2_^2,dz2^x'^)- 



III. COMPASS MODEL 



A simple-looking spin model with similar properties (sometimes called the "compass" model) was introduced in 
1982 in Ref. Q: 

i 

It is very similar to the model Eq. ©, but while the three components of the pseudospin, and (see Eq.©, 

obey the condition + + = , the spin projections S"^, j — x,y, z are independent operators. 

Consider the model Eq.® with J < (ferromagnetic interactions), in which case the angular dependence of the spin 
exchange is similar to that of the conventional dipole-dipole interaction, and indeed, its properties can be qualitatively 
understood by considering the lattice of "compasses" (or magnetic needles). The complications arising in this model 
are immediately clear from this analogy. If, instead of the square lattice, we consider a one-dimensional raw of such 
compasses in the z direction, it is clear that the magnetic needles would order "head to tail" parallel to the z axis. 
However, in a row along the x axis, spins would prefer to order in the x direction. In a square and in a cubic lattice this 
leads to a strong frustration, which in principle can destroy an order altogether. Such claim was indeed made by Glass 
and Lawson for a similar model with real dipole-dipole (long-range) interactions^^: using the Bogolyubov inequality 
they proved the absence of a long-range magnetic ordering at nonzero temperatures in two and three dimensions. 
Interestingly enough, the situation here is in a sense opposite to that in conventional systems: while in standard 
models with a continuous symmetry, like the Heisenberg model, the order is destroyed by fluctuations in one- and 
two-dimensional cases, but can exist in three-dimensional systems, here the order survives in one dimension, but is 
broken for d = 2 and d = 3. This can be traced back to different symmetry properties of the two types of models: 
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FIG. 3: The overlap between d-orbitals of TM ions in the 180°-exchange. 

in contrast to the Heisenberg model, which has a continuous symmetry and, correspondingly, Goldstone modes, in 
Eqs. (|H1) and © the interaction for each pair has an Ising character. However, while in Heisenberg-like models the 
spin and coordinate spaces are largely independent, in the model Eq.|(nil a rotation in the spin space is inherently 
related to the corresponding rotation in the coordinate space. This seems to "restore" the symmetry of the system: 
if we assume a ferromagnetic ordering for d — 2,5 and calculate the spin- wave excitation spectrum, we find that the 
spectrum is gapless despite the discrete Ising type of interactions for each pair. Moreover, this spectrum turns out to 
be one-dimensional for a square lattice and two-dimensional for a cubic lattice: 

Lola^ J^S^ {I -cos k^) (10) 

and 



,l^ = J^S^ {1- cos k^) {1- cos ky) (11) 



^3d 

so that the corresponding fluctuations diverge at all temperatures for d = 2, and at any T > for d = 3 case. 
This agrees with the fact that on the mean-field level there are many equivalent degenerate ground states in our 
system: apart from ferromagnetic states with an arbitrary orientation of spins S, there also states with ferromagnetic 
chains e.g. along the z axis, with spins either parallel or antiparallel to the z-direction and no correlations between 
orientations of spins in different chains. 

What is the true ground state of this model, is still an open problem. One possibility is that beyond the mean-filed 
approximation the order- from-disorder scenario^Jr^: selects a particular ordered state^^. Another alternative is that 
the spin fluctuations are so strong that they destroy a long-range order, leading to a (spin- or orbital-) liquid state^i 
In any case, we see that properties of the models Eqs. ((HJ and are rather unusual. 



IV. EXCHANGE INTERACTIONS IN THE DEGENERATE HUBBARD MODEL 



We would like now to return to JT systems and discuss the purely electronic contribution to exchange interactions 
in these materials. The calculation of the exchange Hamiltonian, resulting from the virtual hopping of electrons on 
neighboring sites, is usually done using the Hubbard model describing electrons on degenerate d-levelsi^ 

H = - ^ ttfdlacrd]f3a + ^ Uafinia^rLip^ ^ S,;i • 8^2 (12) 

ijaf3(7 i.Oi^fi i 

Here i,j are site indices, a denotes spin projection, and a,/3 = 1,2 are orbital indices. The second term in Ea. l|12|) 
describes the on-site Coulomb repulsion between electrons, which depends on the orbital occupation, and the last 
term is the Hund's rule exchange interaction between electrons occupying two different orbitals on the same site. For 
an isolated ion, we have Un — U22 = '^Jh- However, in crystals the direct Coulomb interactions C/q/3 are modified 
by the screening, while the Hund's rule exchange coupling Jh remains unscreened. Therefore, Uafj and Jh should be 
considered as independent parameters. 

For one electron or one hole per TM site, and in the strong coupling limit U 3> Jh ^ W, where W is the electron 
band width, an effective Hamiltonian describing the coupled spin and orbital degrees of freedom can be obtained using 
the perturbation theory in t/ {U, Jh)^- The hopping matrix elements t"j^ are spin-independent (spin is conserved 
in the process of the virtual electron hopping). As a result, the spin part of the exchange operator is rotationally 
invariant and has the standard Heisenberg form Si ■ Sj. In contrast, the hopping amplitudes depend on the type of 
the initial and final orbital, as well as on the spatial orientation of the pair of TM ions (see e.g. FiglSJ, as a result of 
which the orbital part of the exchange Hamiltonian is strongly anisotropic. 
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For the pair of TM sites j and j + z, the hopping only occurs between the two (i3^2_^2 orbitals (see Fig|2K): the 
hopping between the d^z2_r2 and dx2-y2 orbitals is zero by symmetry (see Fig|3lD), and the overlap between the two 
d,j.2_y2 orbitals in the z direction is very small (see FigI3t). Therefore, the usual type of exchange for this pair only 
occurs when the d:}^^^^^ orbitals (T^ ^5) occupied on both sites. By the well-known Goodenough-Kanamori- 
Anderson (GKA) rules, such an exchange is antiferromagnetic, favoring the spin-singlet state of the two electrons 
(holes). This considerations allow us to immediately write an effective exchange Hamiltonian for this pair: 

where t is the hopping amplitude between the two d^z2_^2 orbitals. The exchange Hamiltonian is proportional to the 
product of the projection operators on the orbital states {Sz"^ — r^) on the two TM sites and the projection operator 
on the spin singlet state with Sj ■ Sj^z = The ground state of such a pair is a spin singlet with an orbital 

configuration shown in Fig|21i, in contrast to the (spin-triplet) state shown in Fig|21 stabilized by JT stabilized by 
the electron-lattice interaction. Due to the electron-hole symmetry of the orbitally degenerate Hubbard Hamiltonian 
Ea. (|12|l . the effective Hamiltonians for one hole and one electron per TM site have the same form. There is also 
another mechanism of orbital exchange due to the hopping of electrons on empty orbitals of neighboring sites and the 
Hund's rule ferromagnetic spin coupling^, but it is usually weaker than the antiferromagnetic interaction Ea. (|13|l and 
we shall not discuss it here. 

As in the derivation of Eqs. ^ and the exchange interaction for pairs of TM ions along the x and y directions 
can be obtained from Ea. (|13|l by the corresponding rotation in the orbital space, i.e., by substituting the operators 
in that equation by, respectively, the and 1^ operators. The resulting exchange Hamiltonian for the cubic lattice 
is sometimes called the Kugel-Khomskii (KK) model. The coupled orbital and spin exchange, described by Ea. l|13() 
and the corresponding expressions for Hj^^^ and Hj^j-^-y, has the same spatial anisotropy as the lattice-mediated 
orbital interactions Eq.®, and we face again the problem of frustrated orbital ordering. It was even suggested that 
the coupling between orbitals and spins in Ea. H13|) makes it only worse: it gives anomalously soft mixed spin-orbital 
excitations, resulting in a spin-orbital liquid ground state^Sii^. Other authors, however, argued that this coupling 
actually helps to lift the frustration, resulting in orbitally and magnetically ordered states, in which spins form 
antiferromagnetic chains (as in the cubic perovskite KCuFs), stabilized by large quantum spin fluctuationsif . It 
seems that the only way to clarify this controversy is to perform a direct numerical simulation of the model Ea. (|13|l . 
but the latter is complicated by the sign problem (which again is a consequence of frustration). 



V. 90°-EXCHANGE 



So far, we only discussed a simple cubic lattice of TM ions and showed that even in that case orbital interactions 
are frustrated. There are also many Jahn- Teller materials with lattices that cannot be naturally subdivided into two 
sublattices of TM ions, e.g. a two-dimensional triangular and Kagome lattice, or a pyrochlore lattice of corner-shared 
tetrahedra, e.g. the B-sublattice of spinels. An antiferromagnetic spin ordering in materials with such lattices is 
known to be strongly frustratedi^ and it is interesting to see, if such 'geometric' effect can enhance the frustration of 
orbital interactions. 

The typical example of such a material is LiNi02, in which the low-spin Ni'^+(t2geg) ions with spin S — ^ and 
doubly degenerate orbitals (T = i), occupy two-dimensional triangular [111] layers (see Fig^. While in perovskites 
the angle between metal-oxygen bonds connecting neighboring TM ions is 180°, in LiNi02 this angle is 90°. This 
difference has important consequences for the spin and orbital exchangeiS^. 

The actual exchange in TM oxides occurs not due to a direct d — cZ-hopping, as described by Eq. p2|l , but involves 
the hopping from oxygens to TM sites and intermediate states with holes on oxygen 2p-orbitals. Thus we have to 
start from the d-p model: 

^dp = y ] ^pPima-Pirna- + ^ ] ^dd^j^^jdjacr ~ ^ ] t-ij {p\ma'^jOL(y + d'^j^^Pima^ 
ima jma rnaija 

^ ^ U^Yurt'^^ima^im' 17 ~^ ^ al3^jaa^j (5(7 
imm' JCK/S 

^ ^ ^im ' Sim' Jh ^ ^ ^ja ' ^j(3 (14) 

Here m = x,y, z and a are the orbital indices of, respectively, oxygen 2p-orbitals and TM d-orbitals, and the Hamil- 
tonian of the model includes the direct (Hubbard) and the exchange (Hund's rule) Coulomb interactions on both the 
oxygen and TM ions. 




FIG. 4; The crystal structure of LiNi02. The sublattice of Ni ions consists of triangular [111] layers (shaded). 

In the Mott-Hubbard (MH) region of the general Zaanen-Sawatzky- Allen diagramSli (see also^i), i.e. when the 
charge transfer gap A 3> C/ >• (A is ep — Ed for one hole and Sd + Udd ~ £p for one electron per TM site, and W 
is the band width), the model Ea. (|14|l can be reduced to the effective Hubbard model for d-electrons Ea. (|12|) with 

tdd — On the other hand, in the charge transfer (CT) insulator regime, [/^ ^ A ^ W, one finds new orbital 
exchange terms that turn out to be more important than those obtained within the Hubbard model Ea. H12|l (see the 
discussion below). 

In any case, exchange interactions for the 90°-angle between the metal-oxygen bonds (see FigEJ are very different 
from Ea. (|13|) for the 180°-angle. In the 90°-geometry the orbitals of one TM (e.g. site i in FiglSJ overlap with the 
Pz-orbital of the oxygen 01 (see FiglS)), whereas the c?-orbitals of the ion j overlap with the orthogonal pj;-orbital of 
the same oxygen. Due to the orthogonality of these two oxygen orbitals, it is impossible to transfer a d-electron from 
one TM ion to another via oxygens, and this is the principal difference from the 180°-exchange. The only processes 
allowed are the transfers of oxygen electrons to the site i (from the p^-orbital) followed by the transfer of another 
electron from the Pa;-orbital of this oxygen to the the TM site j. If we first ignore the Hund's rule coupling on 
oxygens, these processes are independent of spins of d-electrons, i.e., they do not result in a spin exchange. However, 
they do depend on the orbital occupation of the two TM ions, resulting in an orbital exchange. Therefore, for this 
geometry, the spin and orbital interactions are decoupled, the orbital exchange being much stronger than the spin 
one (the latter oc Ji//A <C 1). On the other hand, in the KK model Ea. H13|) the strongly coupled spin and orbital 
exchanges are of the same order. Thus, the KK model, obtained for the specific case of the MH insulators with the 
180° metal-oxygen-metal bonds, should not be used uncritically for other situations. 

These general considerations are confirmed by detailed calculationsiS. For instance, for the plaquette formed by 
two oxygens and two TM ions on the sites i and j in the zx plane (see FiglSJ the orbital exchange interaction is given 

by 



Ht — Jt 



n 



(15) 



where 



A2(2A + Up 



A3 



'^'^pd^P 



A3(2A + C/p)' 



(16) 



and the coupled spin-dependent exchange is 



Hts 



231 



(2A + Up) 



(S1S2) 



(17) 



Several points have to be stressed here. First of all, we indeed see that the orbital exchange is stronger than the 
spin one (usually ~ 0.6 — 0.8eV and both Up and A are ~ 4 — 6eV, so that Jf^ < Up, A). The spin exchange 
only appears due to the Hund's rule coupling (in this case on oxygen) and is always ferromagnetic, independent of 
the orbital occupation. Therefore, the spin ordering on this triangular lattice is not frustrated. This, in particular, 
invalidates the attempts to describe LiNi02 as a spin liquid of the RVB typeSSi. 

The second point is the presence of two terms with opposite signs in the expression for orbital exchange Ea. (|15|l . 
The origin of these terms is rather interesting. The negative term resembles those obtained in the Hubbard model: it 
describes the virtual process with two holes on one oxygen in the intermediate state, which is why it contains 2A -I- Up 
in the denominator. The positive term, however, is of a different nature - it remains nonzero even when C/p — > 00. 



FIG. 5: The overlap between the d-orbitals of the TM sites i and j and with two different p-orbitals of oxygen Oi for the 
90°-exchange. 



Its origin can be understood as follows: the virtual hopping of d-electrons from a TM site to 6 surrounding oxygens 
results in an energy decrease. However, each oxygen belongs to two different octahedra. For infinite Up this leads to 
blocking of some of these charge fluctuations: when e.g. the oxygen Oi in Fig|Slis occupied by the hole from the site 
i, it blocks the virtual hole transfer from the site j. This blocking reduces the total energy gain due to the charge 
fluctuations and gives rise to an effective orbital exchange interaction, since an amount by which the energy is reduced 
depends on orbital occupations on the sites i and j. This mechanism of exchange is similar to the Casimir interaction 
between two metallic plates in vacuum. 

While the Hubbard model gives an effective interaction, which respects the electron-hole symmetry, in the d-p 
model this symmetry is lost. In particular, for one hole per TM site the effective orbital interaction the form 



{hole) 



2 ' \2 ^ 



2 ' 



2 



(18) 



(cf . Eq. (|15|l for one electron per site) . In this case the effect of blocking is minimized when the d-holc orbital on at least 
one site is not directed towards the common oxygen, e.g. (i322_r2-orbital on the site j in FiglSja and (ia,2_j^2 -orbital on 
the site j + z. Therefore, the orbital interaction Ea. H18|) favors hole orbitals directed away from the common oxygens, 
in contrast to the situation in MH insulators (see Ea. (|13|l ). for which the energetically most favorable occupation is 
shown in Fig|2^. 

If we now consider the exchange Hamiltonian for the whole triangular lattice, we find out that here too the 
interactions in different directions compete. The outcome is again that in the mean-field approximation the ground 
state is strongly degenerate. The orbital exchange on a triangular lattice is strongly frustrated and has a large number 
of mean-field ground states, including ferro-orbital states, in which (Tj) — Tm on all lattice sites, where m is an 
arbitrary unit vector in the (T^, T^)-plane. Furthermore, there are disordered ground states, which can be obtained 
from the ferro-orbital states by, e.g., inverting the sign of (TJ) on an arbitrarily selected set of lines parallel to the 
lines with the /^/J' type of exchange (see Fig.E)). Such states are ordered along the lines, but there are no long-range 
correlations between the x-projections of pseudospins in the transverse direction (as in the compass model discussed 
in Sec. Illlj) . By circular permutations of the x, y, and z indices one can obtain similar states, which are only ordered 
along the two other sets of lines in the triangular lattice. 

The spectrum of orbital excitations for, e.g., (T) = Tz, calculated in the (pseudo)spin-wave approximation for 
T » 1, is given by 



UJa 



3V2TJj 



sm ■ 



(19) 



As in the case of the compass model discussed in Sec. IIIII the excitation spectrum is gapless and it is independent of 
one of the projection of momentum (q^)- This one-dimensional spectrum leads to divergent fluctuations, to remove 
which, one has to go beyond the (pseudo)spin-wave approximation and take into account interactions between the 
orbital excitations. One obtains then that the order-from-disorder mechanism chooses 6 particular ferroorbital ground 
states, e.g. jSz^ - r^) (all (Tf) = +^)^ and that the spect rum of elementary excitations acquires the gap cx JtVt, 
makes the orbital fluctuations finite. 

We note that the 'geometric' frustration of the triangular lattice plays an important role in this story. One can 
show2^ that, in contrast to the "compass" model, an ordering of Cg orbitals on a square lattice is not frustrated. Thus 
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FIG. 6: A disordered mean-field ground state, in which the pseudospins form hnes parallel to the unit vector exy, such that 
(Tj^) is the same on all lattice sites, while the sign of (T/) varies arbitrarily from line to line. 

in two dimensions, both the anisotropy of the orbital exchange and the special geometry of the triangular lattice are 
necessary for frustration. In a similar way one can treat other types of 'frustrated' lattices. For instance, effective 
interactions Eas. (|15|) and H17|) can be used to describe exchange on the pyrochlore lattice of TM ions, since the 
local geometry of metal-oxygen- metal bonds is the same as in LiNi02. Also the orbital excitations for the d^z^-r^ 
ferro-orbital state on the pyrochlore lattice only propagate along two disconnected sets of bonds with the J^-type 
of exchange, similarly to the case of the triangular lattice. 

Applying the d-p model to systems with the 180°-exchange path, one can show that the 'orbital Casimir' interactions 
also exist in these systems^^. They have the form similar to the phonon-mediated interactions Eq.(|SJl and modify 
the standard KK description. In charge transfer insulators these terms dominate the exchange, resulting in the 
decoupling of orbital and spin orderings and forcing orbitals to order at higher temperatures than spins, provided 
that the frustration of the orbital interactions is somehow lifted. 



VI. SUMMARY 

Summarizing, we have discussed orbital interactions and frustration of orbital ordering in different physical situ- 
ations. We showed that due to the directional nature of the hopping amplitudes the orbital exchange is frustrated 
even for simple (perovskite) lattices. Frustration is therefore a generic property of orbital degrees of freedom. If in 
addition we are dealing with the systems with the 'geometric' frustration, the ground state degeneracy becomes even 
stronger. The resulting type of the orbital, and spin, structure of the ground state can then be determined either by 
small corrections to the exchange interactions lifting this degeneracy, e.g. the Hund's rule coupling or small "external" 
splitting of degenerate e^-levels due to deviations from the perfect cubic symmetry. Another possibility is an ordering 
due to the order-from-disorder mechanism, discussed above. We note, however, that this semiclassical mechanism, 
strictly speaking, works for large values of pseudospin T 3> 1, when quantum fluctuations are relatively small. Since 
in reality T = ^, strong quantum fluctuations may still result in a disordered ground state. Such a quantum orbital 
liquid remains an intriguing possibility. 

The situation with other types of orbitally-degenerate systems, e.g. materials with three-fold degenerate t2g-orbitals, 
may have its own peculiaritics^t. First of all, due to a larger degeneracy and a weaker coupling to the lattice, one 
may expect stronger frustration and stronger quantum effects, which could result in an orbital-liquid ground stated. 
On the other hand, as we already mentioned in the Introduction, in t2g-systems the spin-orbital interaction may 
play an important role, which may remove the orbital degeneracy^^. Furthermore, in frustrated lattices (of the type 
considered in Sec.0) of i2g-ions, one has to take into account not only the exchange between the ions via oxygens, but 
also direct the t2g — t2g overlaps across the TM-oxygen plaquettes (see FiglSJ. What would be the resulting orbital 
structure with all these effects taken into account, is still an open problem. We see, that the orbital physics, especially 
in materials with complex lattices, is very rich and can still produce many surprises. 

This work is supported by the MSC'"'* program. The financial support of the British Council and the Trinity 
College, Cambridge is gratefully acknowledged. 
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